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Abstract 

We consider globally invertible and piecewise contracting maps in higher dimensions and we perturb 
them with a particular kind of noise introduced by Lasota and Mackey. We got random transformations 
which are given by a stationary process: in this framework we develop an extreme value theory for a 
few classes of observables and we show how to get the (usual) limiting distributions together with an 
extremal index depending on the strength of the noise. 


1 Introduction 

In a previous note |7] we announced a few results on the application of Extreme Value Theory (EVT) to one¬ 
dimensional piecewise contracting maps (PCM) whenever the latter are perturbed with a particular kind of 
noise. This noise was introduced in the book by Lasota and Mackey m and it was called randomly applied 
stochastic perturbation (RASP). It turns out that it is particularly useful to study PCM especially because 
it allows very efficient representations of the transfer operator (Perron-Frobenius) and of the stationary 
measures. We now generalize to higher dimensional PCM our previous results by proving the existence of 
extreme value laws (EVL) for some kinds of observables associated to rare visits in the neighborhood of a 
given point. As a matter of fact, we first consider globally invertible maps defined on compact subsets of 
K” and which are diffeomorphisms but on some singular (Borel) sets of zero Lebesgue measure, so they are 
not necessarily contractions. The RASP belongs to the class of random transformations (see next section), 
which naturally define a Markov chain with associated a unique absolutely continuous stationary measure 
(see Th. 10.4.2 in [H]). By taking the direct product of this stationary measure with the probability 
distribution of the concatenated maps, we get the (annealed) probability P which will govern the statistical 
properties of our system. In particular the stochastic process given by the observation along a random orbit 
will be stationary according to P. We will first prove a general result about decay of correlations which 
will not require strong assumptions on the functional spaces of the observables (Proposition 12.II) . We will 
therefore use this result to recover some sort of asymptotic independence in order to prove convergence of 
the distribution of the maxima toward the Gumbel’s law. This will be accomplished in two other steps: 
first we will get an explicit expression of the affine coefficients defining the sequence of levels needed to 
rescale the distribution of the maximum to avoid degenerate limits (Propositions 14.114.214.31) ; then we will 
control the short returns in Proposition 14.51 Whenever we will consider rare visits in the neighborhood of 
a periodic cycle, we will prove the existence of an extremal index less than 1, and in particular equal to the 
magnitude of the noise. 

We would like to stress that we already investigated the EVT for random transformations obtained by 
perturbing with additive noise; we applied it to uniformly expanding maps [T] and even to rotations 0. 
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In both cases we got the existence of an EVL around any point and with extremal index equal to 1. It 
seems that RASP is much more suitable for piecewise contractions and does not smoothen out completely 
the periodicity features of the unperturbed map, as this is reflected in the persistency of an extremal index 
less than one. A few numerical simulations seem to suggest however the existence of an EVL even when 
we perturb PCM with additive noise [7]. 

We now briefly recall the main statements of the EVT. Let us suppose that (Vn)neN is a sequence of real¬ 
valued random variables defined on the probability space (d>, P). We will be interested in the distribution of 
the maximum Mn := maxlLoi Lij ■ • ■ i Vi-i} when n —>■ oo. It is well known that the limiting distribution is 
degenerate unless one proceed to a suitable re-scaling of the levels of exceedances. The precise formulation 
is the following: we have an Extreme Value Law for {Mn)n&i if there is a non-degenerate distribution 
function H :]&—>■ [0,1] with iL(0) = 0 and, for every r > 0, there exists a sequence of levels (u„(r))„gN 
such that 

lim nP(yo > Un) —>■ r, (I) 

n—¥oo 

and for which the following holds: 


lim P(Mji < Un) —>• 1 — H{t). 

n—>-oo 


The motivation for using a normalizing sequence satisfying m comes from the case when 

(Wi)neN s^re independent and identically distributed. In this i.i.d. setting, it is clear that P(M„ ^ u) = 
(E'(u))", being F(u) the cumulative distribution function for the variable u. Hence, condition ([1]) implies 
that 

P(M„ ^ u„) = (1 - P(Vo > u„))" ^ ^ e-", 

as n — >■ cxD. Note that in this case H{t) = 1 — e~'^ is the standard exponential distribution function. Let 
us now choose the sequence = Un{y) as the one parameter family = ?//a„ -I- b^, where y G M. and 
a„ > 0, for all n G N. Whenever the variables V are i.i.d. and for some constants a„ > 0, G K., 
we have P(a„(M„ — bn) ^ y) G{y), where the convergence occurs at continuity points of G, and G is 
non-degenerate, then Gn will converge to one of the three EVLs: Gunibel, Frechet or Weibull. The law 
obtained depends on the the common distribution of the random variables. 

When Yo,Yi,Y 2 , ■. ■ are not independent, the standard exponential law still applies under some con¬ 
ditions on the dependence structure. These conditions will be stated in detail later and they are usually 
designated with D' and D 2 ] when they hold for (yn)neN then there exists an extreme value law for and 
H{t) = 1 — e~'^ , see Theorem 1 in [TT]. We want to stress that these two conditions alone do not imply 
the existence of an extreme value law; they require, even to be checked, that the limit ([T|) holds. It turns 
out that for the kind of observables we are going to introduce, and which are related to the local properties 
of the invariant measure, the limit o is difficult to prove when the invariant measure is not absolutely 
continuous, since one needs the exact asymptotic behavior of that measure on small balls. Instead it turns 
out that whenever the systems is randomly perturbed and the invariant measure is absolutely continuous, 
the existence of the limit (ED can be insured with general arguments and in the majority of the cases (see 
below for a precise meaning), it can be expressed by a closed formula in terms of the strength of the noise. 
Proposition 14.21 and 14.31 below. Moreover that formula could be used in a reversed way: since the sequence 
is now uniquely determined for any n, a numerical sampling for (un)neN which provides conver¬ 
gence to the extreme value law, will bring information on the local geometrical properties of the stationary 
measure: this approach was successfully used, for instance, in Bm- 


2 Random dynamical systems: RASP perturbation 

Let us consider a sequence of i.i.d. random variables (ILfc)fceN with values {ujk)kGN in n space Hg and with 
common probability distribution 9e ■ Let X C be a compact set equipped with the Lebesgue measure m 
defined on the Borel cr-algebra, and {fu)bjens ^ family of measurable transformations such that '■ X ^ X 
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for all CJ € fijil. Given a point x G X and a realization w = {uji,uj 2 , ■ ■ ■) € fig of the stochastic process 
(Wfe)fegN, we define the random orbit of x as the sequence (/”(a;))„gN, where 

f^{x) = X and f^{x) = /„„ o o ■ ■ ■ o f^^{x) Vn ^ 1. 

The transformations will be considered as stochastic perturbations of a deterministic map /, in the 
sense that they will be taken in a suitable neighborhood of / whose size will be determined by the value of 
e, see below. We could therefore define a Markov process on X with transition function 

Le{x,A)= [ lA{fu,{x))dee{uj), (2) 

J Qg 

where A € X is a. measurable set, x € X and 1a is the indicator function of the set A. A probability 
measures is called stationary if for any measurable set A we have: 

= / Ls{x,A)diXe{x). 

Jx 

We call it an absolutely continuous stationary measure, if it has a density with respect to the Lebesgue 
measure. 


We are now ready to introduce the randomly applied stochastic perturbations. We first have a Borel 
measurable map / acting on the compact set X; then we operate an aleatory reset of the state of the 
system at each failure of a Bernoulli random variable: if (x„)„gN denotes the successive states of such a 
random dynamical systems, then at each time n S N we have Xn+i = f{xn) with probability (1 — e) and 
Xn+i = with probability e, where is the realization of a random variable with value in X. This kind 
of perturbation corresponds to the family of random transformations defined by 

fu>ix) =r)f{x) + {l-r])i VxeX, (3) 

where w = ( 77 , is a random vector with value in fig = {0,1} x X. The two components 77 and ^ of w are 
independent and 77 is a Bernoulli variable with the probability of being 0 equal to e, while ^ is a random 
variable that we will suppose Lebesgue-uniformly distributed on A0. The joint distribution 9^ of these two 
components is the product of the Bernoulli measure with weights (e, 1 — e) and the uniform measure on X. 

In order to obtain the stationary measure let us introduce the random Koopman operator [4 : 
L°o ^ defined for all (jx e by 


Ue(j){x) := J (j){fa;{x))d9e. 

Now, if we take two observables </> G and ip G L^, A easy to check that 

J Ue(j){x)ip{x)dx = J J (j){fui{x))lp{x)ddedx -- 


e 


4){x)ip{y)dxdy + (! — £:) 


(p{f{x))ip{x)dx 


4){x)Pe1p{x)dx, 


where is the adjoint operator of Ue, that is the random transfer operator. If we denote P the transfer 
operator associated to / and ip = J ip{y)dy, then we have 


Peip{x) = (1 — e)Pip{x) + eip. (4) 

'^In the following when we will refer to a dynamical system {X, f, p.) we will mean that / is defined on X and preserves the 
Borel probability measure /i; if we will write {X, /), this will simply correspond to the action of / on X. 

^In the book of Lasota Mackay m the variable ^ is more generally distributed with a density with respect to the Lebesgue 
measure; this will not change the results of our paper. 

^From now on and will be referred to the Lebesgue measure m and the integral with respect to the latter will be 
denote as J(-) dx. 
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The stationary measure /ig verifies J (j){x)d^s = J Ue(j){x)d^s and in our case is given by /Xg = /igW where 
/ig G is a density such that /ig = Pg/ig. Such a density exists and is given by |16] : 

CO 

h, = eY,{l-efPH. (5) 

k=0 


In the random setting we are developing the decay of correlations will be of annealed type and it will 
be formulated in terms of the iterates of the random Koopman operator, in particular we have: 

Proposition 2.1. If (j) £ and ■tp £ D then 


J J (I^{f^{x))'fp{x)dfied0^ - J (j){x)dne J i’ix)dn. 


^ 2(1 — £)"'||^||loo Il'i/j/igll 


LI- 


( 6 ) 


Proof. To use the duality of the Koopman operator with the transfer operator P, let us introduce the 
following quantity: 


Corm{4>,ip,n) := 


Ufp{(j){x))ip{x)dx — J (j){x)d^e J 'tp{x)da 


(f>{x)Pf‘'tp{x)dx — J (j){x)d^e J 'ip{x)dx 

(7) 


By recurrence, it is easy to obtain the following formula for the iterates of P^ : 

n —1 

P"i/; = (1 - £)"P"V' +efp'^{l-e)'^p'^l. 


/c =0 


By replacing this expression into dZD and writing there d/ig = h^dx, with given by ([S]), we get 


Corm{4',i>,n) = 


n n OO 

(1 - erj 0(x)p-^(.)dx - .yy ^(x) i:(l - x)*P‘l(x)<i. 

k—n 


< 2{l- eYWPWl’^W'PWli, 


k—n 


( 8 ) 


by Holder inequality and the fact that the infinite sum is of positive terms. Therefore, it ip £ n 
then the correlations ([5]) are equal to 


Corm{(p,iphs,n) ^ 2(1 - £)”||(/)||L=o||^/ig||ii. 


□ 

Remark 2.2. ITe will argue later on that in the EVT the preceding decay of correlations result plays a 
fundamental role whenever ip is a specific characteristic function: in this case it will be sufficient to have 
the density he in ffi. For observables ip which are solely in ffi we will need a density in L°°. Note that 
the exponential decay of the correlations does not depend on the map f, which is a remarkable property of 

RASP. 

We are now ready to introduce the extreme value statistics; at this regard we define a stochastic process 
(Ki)„gN, n G N, by composing a given observable <p : X ^ M. a.s: Yn = cpo /”, with 

=T]nf{x) + {l-r]n)f,n, V X G X, 

where a;„ = {rjnffin) £ {0,1} x X is random variable which distribution is the product of the Bernoulli 
measure with weights {e, 1 — e) and the uniform measure on X. 

In this case (Kn)ngN will be a stationary process if we consider the probability P = pg x 9f, which 
corresponds to the annealed situation where we average over the initial condition and over the realization 
of the noise. 
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We will in particular choose the observable ^(a:) = —\ogd{x,z), where denotes (some) distance 

on X and z a point oi X. As we anticipated in the Introduction, such an observable is related to recurrence 
in small sets, since the distribution of the maximum of the random observable (j)of^^k = 0,...,n — 1 up 
to the level Un, coincides with the distribution of the first entrance of the random orbit into the ball (in 
the metric d), B{z,e~‘^’') centered at z and of radius e““" 

We now state the two conditions which insure weak dependence of the process and which allow to get 
the limiting distribution of the maxima. 

Condition[Z32(u„)] We say that D 2 {un) holds for the sequence Iq, W, • ■ • if for all £, t and n. 


|P (Fo > Mn n max{Ft,..., Ft+^_i ^ «„}) - P(Fo > u„)P(M^ < Mn)| ^ l{n, t) 


where 7 (n, t) is decreasing in t for each n, and n'^(n, f„) —> 0 when n —>■ oo for some sequence tn = o{n). 
Now, let {kn)nGN be a sequence of integers such that 

^ oo and = o{n). (9) 

Condition[ZI'(u„)] We say that D'{un) holds for the sequence Yo,Yi,Y2,... if there exists a sequence 
(fcn)raGN Satisfying dH]) and such that 


ln/k„\ 

lim n P(Fo > u„,F > u„) = 0. (10) 

i=i 

As we said in the Introduction when these two conditions hold for (Fn)nGN then there exists an extreme 
value law for and H{t) = 1 — e~'^ , but provided that the limit ([1]) is true. 


3 Stationary measure for injective piecewise diffeomorphisms with 
RASP 

We now introduce the dynamical system upon which we will perform the randomly applied stochastic 
perturbation. 

Definition 3.1. The map. Let us X be a compact subset of which is the closure of its interior 
X = int{X), equipped with the metric d and with the normalized Lebesgue measure m; suppose is a 

collection of N disjoint open subsets of X such that X = ur= Xi and m(A) = 0, where A := X \ Xi 
is the singular set. We will consider non singular maps f : X ^ X which are injective and such that /|x\A 
is a -diffeomorphism. We say that f is continuous in a point w of the boundary of X if there is an open 
ball of radius k, B{w, k) such that f is continuous on B{w, k) H int{X) and it can be extended continuously 
on B{w, k). 

Remark 3.2. In the examples we will treat, the set A will be composed by points where f is discontinuous 
and by the points of the boundary of X. We would like to stress that a point of that boundary is not 
necessarily a discontinuity point of f. 

The following proposition give the expression of the stationary density. 

Proposition 3.3. Let {AfcjfegN be the sequence of sets defined by Aq := X and A^+i := /(A^ \ A) for all 
k ^ 1 . Denote Jk{x) := OiLi I bet(/'(/“*(a;))|“^ for all k ^ 1 and Jo{x) := l(a;). Then, 

CO 

hs{x) = - £)'^Jk{x)lAk(.x) \fxeX. (11) 
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Note that Jk{x) is well defined for every x € X such that x ^ /*(A) for alH G {1,..., fc}. Among other 
basics properties of the sets which will be helpful all along this paper, the following lemma ensures that 
Jk{x) (and therefore he{x)) is well defined for all x G Ak- 

Lemma 3.4. The sequence of sets {Afej^gN has the following properties: 

1. Ak is open for all k ^ 1 

2. Afe_|_i C Afc for all A: G N, 

S. Afc n /P(A) = 0, for all p G {I,..., k}. 

Proof. 1) The set Ai is open since X \ A is open and / is a diffeomorphism in X \ A. Using the same 
arguments, we obtain by induction that Afc is open for any k ^ 1. 2) We have A 2 = /(Ai \ A) C f{X\A) = 
Ai and if we suppose Afc C Afc_i, then Afc+i = /(Afc \ A) C f{Ak-i \ A) = Afc. 3) By injectivity we have 
f{X \ A) n /(A) = 0, i.e. Ai n /(A) = 0. Suppose now that Afc n /^(A) = 0 and that there exists 
X G Afc+i n /^+^(A). Then there exist y G Ak\ A and y' G f^(A) such that f(y) = f{y') = x. As / is 
injective, we must have y = y', which contradicts the recurrence hypothesis. It follows that Afcn/^(A) = 0 
for all k ^ 1 and applying 1) we obtain the desired result. □ 

Proof. (Proposition [331) • We compute the transfer operator P of / in order to use the formula ([S]). Since 
/ is a diffeomorphism on the full Lebesgue measure open set X \ A, its associated transfer operator acts 
on a function ip G L^{X) as 

Pip{x) = lf(x\A){x)ip{f~'^{x))\det{f{f~'^{x))\~'^ = lAi(a;)^/(/"^(a;))Ji(a;) Vx G A. 

Since the set /(A) has zero Lebesgue measure, we can suppose without loss of generality that Pip(x) = 0 if 
X G /(A). Then, Pip{x) is well defined for all x G X, since f{X \ A) = Ai does not contain point of /(A). 
Now, we show by induction that for the uniform density 1 = lx we have 

P^l{x) = lAk{x)Jk{x) V/c > 1. 

It is easy to check that Pl(a;) = 1 ai (a;) Ji(a:) and that 

P'=+H(a:) = lAAx)lf(Ak)ix)Jk+i{x), 

if we suppose the property true at rank k. Now, we have 

Ai n /(Afc) = (Ai n Afc+i) u (Ai n /(Afc n A)) = Afc+i u (Ai n /(Afc n A)). 

Since Ai n /(Afc n A) c Ai n /(A) = 0, we have Ai n /(Afc) = Afc+i. It follows that P'"'+^l(a;) = 
lAs,^_i(a:) Jfc+i(a:), which ends the induction. Using ([S]) we obtain that 

00 

he{x) = - e)’^Jk{x)lAk{x) VxGA, 

fc =0 

which ends the proof of the proposition. □ 

Now we decompose the space X as the disjoint union A = Ai U /(A) U (A \ /(A)) and we use (fTTll 
to study the density in each of these sets. Since the last two sets do not belong to any Afc with fc 7 ^ 0, 
we obtain that h^ = e on these sets. For a; G Ai, by property 2) of Lemma [3.41 two situations can occur: 
either there exists p ^ 1 such that x G Ap \ Ap+i or a: G A := Afc. So in Ai we have on one hand 

p 

heix) = - e)^Jk{x) Va:GAp\Ap+i, (12) 

k^O 

and the by continuity of the the density is continuous. On the other hand we have 

00 

he{x) = — e)'^Jk{x) x G A. (13) 

fc =0 
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In this last case the density may not be bounded, depending on the behavior of Jk{x). However, if we 
suppose that there exists 0 < A < (1 — e)~^ such that \det{f'{x)\~^ ^ A for all cc € X \ A, then the series 
(I13p converges uniformly in A and is also continuous in A and thus on the whole set Ai (see Proposition 

m. 

In the next sections we decompose our study of the existence of EVL in two parts. The first part 
concerns the point of Ai \ A, where 

A := Pi Afc. (14) 

feeN 

The set A is called the global attractor of / and contains all the point of A. In Ai \ A we can verify the 
conditions D 2 {un) and D'(un) to show the existence of EVL, but in A condition D'{un) is not verihed. 
For this reason in the second part, we introduce other conditions which allow us to show the existence of 
EVL but with an extremal index different from I; this is achieved by giving additional property, namely 
by choosing / as a piecewise contracting map. It has been shown in [5] that for any piecewise contracting 
map defined in a compact space, if the global attractor A does not intersect the set of the discontinuities, 
then it is composed of a finite number of periodic orbits. Moreover, this condition, for injective maps as 
our, is generic in topology [5]. All these properties will help to get the existence of EVL in the set A. 
Let us give a few example to illustrate that 

Example 3.1. As an easy example of an uncountable hyperbolic attractor, let us consider the baker’s map 
f defined on the unit square [0,1] x [0,1] iteratively for n > 1, by 


Vn 


^ f ivn, for yn<a 

\ J^ivn -a)+ jbXn, for yn> a 


with 0 < Xn, < 1, 0 < 7 o < 7 & < ^, a < i. /n this case A is given by the union of the boundary of the 
square and the segment y = a. Moreover Ai is given by the two rectangles [0,7a] x [0,1], [|, ^ +7b] x [0,1]. 
Notice that the derivative f is a diagonal matrix with diagonal elements 7 q (or^b), ond ^ (or y^)- 
can choose these parameters in order to get det(/') < I. 


Typical examples of PCM for which the attractor is generically a finite set of periodic orbits are piecewise 
affine Him. That is, if for each i G {!,...,A}the restriction fi := f\xi of the map f : X ^ X of Definition 
o to a piece Xi has the following form: 


fi{x)=AiX + Ci VxGXi, 
where Ai : K.^ —>■ is a linear contraction and Ci € X. 

Example 3.2. 1-D case The simplest example of a piecewise contracting map of this class is given in 
dimension I in the interval [0,1] by f{x) = ax + c with a,c € (0,1) (here Xi = (0,(1 — c)/a), X^ = 
((1 — c)/a, 1) and A = {0,1, (1 — c)ja\). For almost all the values of the parameters a and c, the global 
attractor A is composed of a unique periodic orbit (which period depends on the specific values of the 
parameters), but for the remaining set of parameters the global attractor is a Cantor set supporting a 
minimal dynamics. We refer to for a detailed description of the asymptotic dynamics of this map. We 
remind that in our paper m we investigated in detail the case with c = 0 and 0 the (unique) fixed point; we 
briefly quote here the simple structure of the density 

he{x) = eX! ^x e{aF,aF~\ p>l 

and the formulae for the coefficients 

Un = F bn with 

an = 1 o,nd 


of the linear scaling. At this regard let r > 0, y = — ln(r) and 

p-i 


bn = log 


l2neY, 


(i-^r 


Vn G N, 




7 




(i) z 0 on the interval but not in the countably many discontinuity points of hg, namely z ^ 

(ii) p'^l such that z £ and n is large enough such that the ball _B(z, e““") C then: 

nF{XQ >Un) =T 

Example 3.3. 2-D case In the unit square X = [0,1]^, another interesting example is given by a map with 
four pieces Xi := (0,Ti) x (0,T2), X2 := (Ti, 1) x (0,T2), X3 := (0,Ti) x (T 2 , 1) and X4 := (Ti, 1) x (T 2 , 1), 
singular set A given by the two segments {x = Ti}^{y = T 2 } and the boundary of the unit square, and 
restricted maps fi{x,y) := a{x,y) + (1 - a)(l,0), f 2 [x,y) := a{x,y) + (1 - a)(l,l), f 3 ix,y) := a{x,y) 
and f 4 (x,y) := a(x,y) + (1 — a)(0,l). Once again, for almost all values of the parameter a £ [0,1), Ti 
and T 2 , the global attractor is composed of a finite number of periodic orbits. However, contrarily to the 
1 dimensional case, where the number of periodic orbits is bounded by the number of pieces N m, the 
number of periodic orbits increases with the parameter a and tends to infinity when a goes to 1. This 
example belongs to a larger class of higher dimension piecewise affine contracting maps which are models 
for genetic regulatory networks and are studied in m- 

Question An interesting question is the check if the stationary measure is stochastically stable, in the 
sense that the measure fie will converge weakly to the SRB measure supported on the global attractor. Let 
us remind that it is not the case for the simple contraction f{x) = ax, x £ [0,1], a < 1, as it was proved in 
Lasota Mackay m- 


4 Convergence to the extreme value law 


We prove in this section the convergence towards an extreme value distribution in the complement of the 
attractor defined in (fT^ . We will work with the process W(x,w) := — log d{f^{x), z). Our first task will 
be to prove the existence of the limit O; this will be accomplished in two manners: we will first give a 
general result showing the existence of that limit for any point of X, but without an explicit expression for 
the scaling coefficients a„ and Such coefficients could instead be computed on the set Ai \ A and even 
on the attractor itself provided the density of the stationary measure will be essentially bounded. The first 
general result uses Theorem 1.7.13 in the book by Leadbetter, Lindgren and Rootzen [15]: according to it, 
a sufficient condition to guarantee the existence of the limit (U]) for 0 < r < 00 , is that Ij 

X —>■ xp, where F is the distribution function of Yq, the term F(x—) in the denominator denotes the left 
limit oi F at a: and finally xp := sup{a:;E(a::) < 1}. For our particular observable To, xp = 00 and the 
probability P reduces to the stationary measure yie which is absolutely continuous and therefore non atomic. 
Therefore F is continuous and the above ratio goes to 1. We have thus proved the following proposition: 

Proposition 4.1. Suppose f X ^ X satisfies Definition \3.1\ Then, for any z £ A defining the observable 
^o(') = ~ z) o.'<xd any 0 < r < 00 , there exists a sequence Un such that 

lim nP(Ao > u„) = r. 

71^00 


As we anticipated above the previous result could be strengthened whenever z lies in the complement 
of the global attractor under the general assumption that h^ is summable. Actually, we will suppose that 
z £ (Aq \ A) \ where 9Afc := Afc \ Afc; this in turn implies that z £ Ap \ Ap_|_i C Ap \ Ap+i for 

some p ^ 10 We can take advantage of the continuity of the density in such points to obtain a possible 
scaling sequence for which the limit (IT|) holds. 

Proposition 4.2. Suppose f : X ^ X satisfies Definition \3.1\ Let p ^ 1 such that z £ Ap \ Ap+i. Let 
T > 0, y = — ln(r) and Un = -^ + bn with On = D and bn such that 

m{B{z, = ins 


"^We notice that in the case of the baker’s transformation of Example 1 3.1 1 the boundaries of the Afc belong to the attractor 
and therefore we can take simply z G (Aq/A). 


j2i^-sfMz) 


k=0 


nhs{z) 


Vn £ N, 





where B{z,r) denotes a ball of center z and radius r. Then, 

lim nP(Xo > Un) = t. 


(15) 


n—?'00 


Proof. For all n S N, we have 

nW’{Xo>Un) = n j = n J lB{z,e-'^n){x)he{x)dx. 

Since (Mn)ragN is an increasing sequence for any r > 0 and 2; is in the open set Ap\ Ap+i, if n is large enough 
C Ap \ Ap+i. Using formula (fT^ and m{B{z,e~'^'')) = Tm{B{z,e~^"-)) we obtain: 

p p 

nTm{B{z,e~^^))e'S^{l—e)’^ inf Jk{x) ^ nF{Xo > Un) ^ nTm{B{z,e~^^))e'S^{l—e)’^ sup Jk{x). 


fc =0 

Since Jk is continuous in B{z, e““”) for all k G {0,... ,p}, we have 

lim inf Jkix) = Jk{z) = lim sup Jk{x). 

n-*-ooxGB(z.e-“n) a:GB(z.e-“" ) 

Using the expression of the sequence (Mn)nGNi we obtain 

lim nP(Alo > = r. 


(16) 


□ 


We note that an inspection of the proof of Proposition 14.21 easily shows that similar results can be 
obtained for the points of A, but supposing a continuous bounded density. The following proposition shows 
that the density has such properties in A under some condition on the derivative of /. 

Proposition 4.3. Let 0 < A < (1 — e)“^ and suppose that \ det(/'(a;)|“^ < A, for all x G X/A, then is 
continuous in A and the same conclusions of Proposition \4-S\ hold. 


Proof. For any k > 0 and a: G AT let gk{x) := e(l — e)^ Jfc(x)lAj,(a:). Then, gk{x) = 0 if a; A^ and 
gk{x) < £:((! — e)A)^ if x G Aj, and therefore Yl^k=o9k{x) converges uniformly on X when p —>■ 00, since 
^^o((l “ ^)''')^ converges. Now, let xq G A, then 


lim he{x) 

X—^Xq 


lim gk{x) = 

cr. — ^cr.r\ 


/c=0 




lim gk{x) = £^^{1 - e)^ lim Jfc(x) = e 'V(l - e:)'" Jfc(xo) 

X—¥Xq X—^Xq Z—/ 

/c=0 fe=0 


he{xo) 


We have used that xq G A^ for all fc > 0, which implies on one hand by openness of the Afc that 
limx^xo^Akix) = 1 and in the other hand that Jk is continuous in xq, since for any x G A^ we have 
/*(x) ^ /(A) for ain < fc and / is in Af \ A. □ 


We are now ready to check conditions 112(u„) and D'{un) which will guarantee the convergence of our 
process toward an extreme value law, the Gumbel’s one in our case. As we said above condition D 2 {un) is 
insured by the decay of correlations result established in Proposition 12.11 see the proof of Th. D in Sect. 
4.1 in [1] for the very general approach. We are thus left with the proof of D'{un). We first need: 

Lemma 4.4. Let n G N and Un = B{z,e~^’') be a ball centered at a point z G X with radius e““". Let 
p ^ 1 be such that z G Ap \ Ap-)-i then, 


min{j-l,p} 


.. ,ujj) : X G UnJiix) eUn) = efieiUn) ^ (l “ / Jk{x)lur^{x)dx Vj > 1, (17) 


for n large enough. In particular, under the same hypothesis we have, 

F{{x,UJi, . . . ,UJj) : X G Un,fi{x) G U„) ^ /ie(C/„)^ Vj ^ 1. 


(18) 
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Proof. Let us denote 


Prj :=P((a;,a;i,...,u;j) : x & Un,fi{x) G t/„) Vj > 1. 
First of all, for all j ^ 1 we can write 


Pr,= 


J J lc/„(a;)l(7„(/i(a;))d6»^(i/re = J lu„{x)Ui{lu„)ix)dne = J Pi{lu„h^){x)lu„ix)da 


Using ([8]), we obtain that for all j ^ 1 

Pr^ = [ PH^u„K)ix)lu„{x)dx 


{1 - ey J P^{lu^he)ix)lu^{x)dx + e^ieiUn)'^{l - ey J P''l{x)lu„ix)dx 

k—0 

(1 -ey f lu^{x)lu^{P {x))he{x)dx efj.e{Un)'^{l -e)^ /* J/c(:r)lAfe (x)la^ (x)dx 
{I - ey f~yUn)) + ^^ie{Un)^{l - / Jk{x)lK^nUyx)dx. 

1 -n 


Now suppose n large enough for Un being a subset of Ap \ Ap+i. Then one can show that fyUn) C 
U'fc=i /^(^) 'J ^p+j (using recurrence on j). Now, in the union [yk=i f^i^) have that Ap n /^(A) = 0 
for any k ^ p. We also have /^(Ap) n /^+*(A) = 0 for any I G N, since Ap n /^(A) = 0 and / is injeetive. 
Together with Ap_|_j C Ap, this implies that Un H A(t/„) = 0, and therefore C/„ H f~yUn) = 0 for all j ^ 1. 
On the other hand, Un Pi A.k = 9 for all k > p. It follows that 


min{j-l,p} 

Po = ^dXUn) Y1 “ y'" / 'h{x)\uyx)dx Mj > 1, 

k=0 

which prove (jl7l) . Now, it easy to show using (fT^ that Pr^ ^ Prp+i = for any j ^ 1. □ 

Proposition 4.5. Letp ^ 1 and z G Ap\Ap_|_i, then the maxima of the proeess Xn{x,w) := — ln(i(/"(a:), z) 
verifies condition D'(un). 


Proof. Suppose z G Ap \ Ap+i for some p ^ 1. Consider the sequence (u„(?/))neN of Proposition 14.21 Let 
{kn)neti be such that —?► oo when n —?> cxd and = o(ri). Supposing n large enough and using (fT51) . we 
deduce that 


LaJ 

n ^ P((a;,wi,... ,a;j) : x G Un,f^{x) G Un) ^ n 




(npe(Un))^ 

kn 


Using \im.np^(Un) = t, we obtain that 


lim P((a;,wi,..., Wj) : x G A o • • • o f^^x) G Un) = 0, 

n—^oo 

i=i 

since np^Un) —t r and A:„ —>■ oo when n oo. In such a case, (A„(x,w))„gN verifies conditions I? 2 (it„) 
and D'(un) for our sequence {un{y))n&ii wich implies that P(M„ ^ Un{y)) converges towards the Gumbel 
lawe-®”". □ 


Collecting all the results of this section we finally proved the following theorem 

Theorem 4.6. Suppose f : X ^ X satisfies Definition [SAi Then, for any z G (Ai \ A) \ U'^ 2 ^^k, the 
sequenee {Mn}„gN of the maxima of the proeess defined for every n G N by Yn(x,ui) := — log (i(/”(x), z) 
admits the Gumbel’s law as extreme values distribution. 
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5 Clustering in the attractor of piecewise contracting maps 


In this section, we complete the results of the results of Theorem 15.11 by studying the case when the 
maximum of the observable is achieved in the set A of a piecewise contracting map. This creates clustering 
of exceedances which are responsible for the appearance of an Extremal Index less than 1. To prove the 
existence of such an EVL, we follow the techniques developed in m that were recently upgraded in la¬ 
in particular, we will use the notation used in the latter. Actually in m the authors introduced two 
conditions on the dependence structure of the stochastic processes in order to obtain convergence to an 
EVL under the presence of clustering. These conditions were called DP{un) and D'^{un) and essentially 
were adaptations of conditions Z? 2 ('a„) and D'{un), specifically designed to cope with clustering. Later, in 
|14| . the authors introduced some general conditions called JJ,{un) and which polished and joined 

all the previous conditions, allowing to address the presence {q ^ 1 ) and absence {q = 0 ) of clustering at 
once. We recall this conditions here. 

Let us begin to introduce a few notations. For q gNq we put 

^ := {^0 > Un, Yi <Un,...,Yg < M„} 


where again Yk,k S Nq is a P-stationary process, that we will take as that studied in the preceding chapters. 
For s,£ G N and an event B, let 

s+i-l 

WsAB)= n (19) 

i—s 


Condition We say that holds for the sequence Toj Ti,... if for every £,t,n € N and q € Nq, 

P n (a^)) - P (a^) P (Wo,i (a)^))) I ^ ^iq,n,t), (20) 


where 7 ( 5 , n,t) is decreasing in t for each q,n and, for every q G No, there exists a sequence (tri)nGN such 
that tn = o(n) and nj{q,n,tn) —>■ 0 when n —> c». 


For some fixed q G Nq, consider the sequence given by condition /l,(ti„) and let (A:„)„gN be 

another sequence of integers such that 


^ 00 and kntn = o(n). ( 21 ) 

Condition (,Z(^(m„)). We say that JXg(un) holds for the sequence Yg, Yi, I 2 , ■ • ■ if there exists a sequence 
{kn)neN satisfying (| 21 |l and such that 

[n/fcrij —1 

Jh^ n ^ P (a(«) n T-^ (A(f))) = 0. (22) 

i=9+i 


Now let 


1 } = lim -dn 

n-s-oo 


lim 

n—>-oo 


P(a1«^) 
nUn) ■ 


(23) 


From [H Corollary 2.4], it follows that if the stochastic process Yq,Yi,. .. satisfies conditions JJ,{un) 
and P^{un) and the limit in (l23l) exists then lim„_>oo ]P(AI„ ^ «„) = where the sequence (M„)„gN is 
such that lim„_,.oo nF{Xo > u„) = r. 


5.1 Analysis in periodic attractors 

Using the result just mentioned we will show that 

Theorem 5.1. Suppose that f : X ^ X satisfies Definition \3.1\. Suppose moreover that f is piecewise 
contracting, that is, there exists a G (0,1) such that for every i G {!,..., N} we have 

difix), fiy)) ^ adix, y) Vx, y G W. 
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Suppose that the set A does not contain any discontinuity point of /, and let z G A. Consider the random 
perturbations defined in ([3]) and let ¥ = , where /ig is a stationary probability measure. Let Yq, Yi, ... 

be the stochastic process given by Yn{x^uf) := — logd(/”(a:), 0{z)), where 0{z) is the orbit of z. Let (un)neN 
be a sequence such that lim„_>.oo n¥{Xo > Un) = r ^ 0. Then we have 

lim ¥{Mn ^ Un) = e~^'^. 


It has been shown in [1] that for any piecewise contracting map defined in a compact space, if the global 
attractor does not intersect the set of the discontinuities, then it is composed of finite number of 

periodic orbits. Therefore, under this hypothesis 0{z) is a periodic orbit of the map. If p is the period of 
this orbit, then {Yq > Un} occurs if and only if x € C/n, where 

p-i 

Un= Us(r(z),e-“"). 
i =0 

Remark 5.2. Note that on the contrary to the usual choice for the observable function, here, we are consid¬ 
ering that the maximum is achieved at every point of the whole periodic orbit, instead of at a single periodic 
point. In an ongoing work concerning EVLs for observables achieving a maximum on multiple correlated 
points, the second-named author with D. Azevedo, A.C.M. Freitas and F.B. Rodrigues are developing a 
thorough study regarding the consequences of having multiple maxima in terms of the extremal behavior of 
the system. The important fact to retain here is that by considering the whole orbit we are performing a 
sort of reduction to the case where the period is 1. 


As in the previous section. Proposition 14.11 guaranties the existence of the a sequence {unjneN such 
that the limit (HD holds. However, as a side comment, we note that we can obtain an estimation of the 
sequence using a similar proof as the one of Proposition 14.21 provided the density is bounded and 

continuous. Since A n A = 0, it is composed of finite set of periodic points that are not in A and using the 
definition of the sets A^, it is easy to show that all the points of A also belongs to A, i.e A = A. Thus, we 
can argue that, under the hypothesis of Theorem 15.11 and of Proposition 14.31 the density has the desired 
properties on A to obtain an estimation of the sequence {unjneN- 

The proof of the theorem consists in showing that Foj • satisfies conditions JXfun) and JXq{un) and 
moreover d given in (1^51) is well defined and equals e. Since 0(z) n A = 0, for n large enough each ball of 
the union defining t/„ is contained in one of the open pieces Xi, and because of the contraction we have 
that f{Un) C Un- Therefore, we will consider q = \ and use An'^ = {(x, w) : x G f/„, /^(x) ^ Un}- 


Verification of 

From the proof of Proposition 12 .1 1 we have: 


Corm{4',tp,n) := 


JJ 'fp{x)(j){f^{x))de^{u)dx - J tlj{x)dx J (l>{x)dne{x) 


<2(l-£r||</.|U||V’lli- (24) 


We will use this information about decay of correlations to verify condition in the case <7 = 1. We 

have 



+ (l-e )II ^u„{x)lu-{f{x))l^^^^^^ii),^{fl'^^{f{x)))d9^{u)dp.e{x) 

= £Re{Un) II 

= efi,{Un) II Iu^{y)(l>{flz^{y))d9l-^{UJ2,... ,ajt)dy, 
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where 


(l){x) = J ... ,uii). 

In the third equality above we used the fact that if f{Un) C Un and consequently lun{x) ■ Iw^ifix)) = 0 
for SiW X & X. We will use this fact again in the following computation: 

= JJ '^uAx)'i-u-{fu;J{x))de^{^)dne{x) = e jj lu„{x)lu-{y)dyd^Xe{x) 

+ (!-£:) J lu^{x)lu^{f{x))dys{x) = eye{Un) J lu-{y)dy. (25) 

Also 

= jj l^^ ^^^(^)^{x)ef{ui)dye{x) = J 4 >{x)dye{.x). 

It follows by ((24l) 



Hence, we may take 


— ^deiUn) 


'^m{y)<i}Uauj iy))dd,i^)dy - / iu-{y)dy / 4’ix)dyeix) 


= £ye{Un)Corm{(j), - 1) ^ 2epe(C/„)(I - ef 

logn, for example, and is easily verified with = 2£fig(Un){l — 


ey. 


Verification of P^{un)- 

For j > I, let Pr^^^ = P({(x, w) : x € A^y fi{x) € We have 

^ JJ ^u„{x) ■ l(7c(/^(x)) • luAfiXx)) ■ lc/c(/^+^(x)) dO^ he{x) dx 
III ■ ^uAfL^iy)) • ^u-ifL{y)) dy de^ h,{x) dx 

+ (1 -e) JJ '^uAx) ■ lu-ifix)) ■ iuMaZ^ifix))) ■ lu-ifLifix))) dey h,ix) dx 
= ^ J 'i-uAx)he{x) dx Jj 1(7= (y) • l(7„(/;^^^(y)) • 1(7= (/^^(y)) dy dO^ 

< sy,{Un) Jj lc/=(y) • luMiZ^iy)) dO^ dy = ey,{Un) Jj l(7=(y) • C/J“^(lc/„)(y) dy 
~ (Ufi) II Pi ^(ic/=)(y) • ic/„(y) rfy, 

where in the third equality we used the fact that if x S t/„ then /(x) ^ C/^ and consequently Xjj^ix) ■ 
1 ( 7 = (/(x)) = 0 for all X & X. 

Now, using ([8]), we define: 

j — 2 

JJ Pr\'^u-)iy) ■ 'i-uAy) dy = {I - £y~J P^~\lu-)iy)'i-uAy)dy + em(u^)^{l - e)J p'=l(y)lc/„(y)dy 

k—0 

J-2 

= (1 - ey-^m{U^ n f-yUr,)) + emm ^(1 - s)'=m(/-'=(C/„)) 

k^O 

Note that to check we need to show that n Prj^^ vanishes. But n Prj^^ ^ 

n£iJLe(Un)Yl}jll 2 ''^ ^ k + Pp Since Imin^ao nfj,e(Un) = T > 0, then to verify JX^Un) we only need to 
check that both IIj vanish as n —>■ oo. 
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Since A n A = 0, there exists /cq := niin{fc € N : n A = 0} and since the sets A^ are open and 
satisfy A^+i C A^ for all /c G N, for a n large enough there is a fc' ^ fco such that Un C A'^,. Therefore, 
we can define dn ■= ko + d'^, where d'^ := max{fc G N : 17„ C A^}. Note that liin„^oo dn = oo, since the 
distance d{Ak,A) —>■ 0 when k oo. Moreover, as for any k ^ ko we have Afc+i = /(Afc), it follows that 
f~^{Un) C f~'^{Ako+d„) = for all k G {0,1,..., d„}. Also, for any k G {0,1,..., d„}, we have 


^0"r“Ti,/ J \ >^0) 

P ^ (T T \\ ^ ( -P^n ^ / A. ^ ^ / 


m{f (t/™)) < ^{Ako))^a '‘/3„ where /3„ := a‘^"m(AfeJ. 

Now we can check (1^^ starting with splitting the sum 

[n/kn]—! ln/kn\—2 [n/fc^J—2 

= X! ^i+i=x!^*+i+ X! 

j — 2 2 = 1 2=1 2 = <in + l 


For the second sum on the right we have: 

Ln/fe„J—2 oo 

( 1 -e)* ^ -e)" 

i—dn-\-l i—dn-\-l 

For the first sum on the right, observe that 

dn dn dr, 


t + 1 


->■ 0 , because lim„_>oo dn = oo. 


^ A+i = ^(1 - eYmiU^ n rYUn)) ^ y](l - eymif-YUn)) ^ 


2=1 


2 = 1 


Now to analyze the asymptotic behavior of /5n X!^=i ~ ^ consider three different cases. We 

assume first that (1 — £)a“^ > 1 . 


/3„ ^(1= /3, 


[(1 — e)a — (1 — e)c 

(1 — e)a“^ — 1 




[(1 — e)a 

(1 — e)a~^ — 1 


^ TO(Afco) ^ lim„_j.oo dn = oo. 

yi — e) — a n-j-oo 


In the case (1 — e)a ^ < 1, then 

dn OO 

/3n y^(l - e)*a“* < l3n y](l - efc 

i=l i=l 

In the case (1 — e)a~^ = 1, then 


0 , since lim„^oo Pn = 0 . 


dn 


Pn ^ (1 ^') Yy ^ ^ dri^ ^ 0 ; since lim ^^—^.oq dji — oo. 

• ^ n—>-oo 


2=1 


Now, we turn to ^ Let pn and dn be defined as above, then 


[n/fcnj—1 [n/fcnj—3 dn 2 [n/fcnj—3 i 

IIj= E II^+2=^em{U^y)Y,i^-e)>^m{rHUu))+ E - e)’^m{r\Un)) 

j — 2 2=0 2 = 0 k — 0 2=rfn + l fc=0 

dn 2 L'^/^nJ OO 

< EeE(^ + E 


i—0 k—0 


^ e^n ( 1 + E E*-^ “ 


2=1 k—0 


2=1 k—0 
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Here we have used (O to obtain 


/ OO « OO 

^u„ix)he{x)dx = e'^il - e)^ / P''l(a;)l£/„(x)dx = e^(l - £:)''m(/“''(C/„)). 
fc =0 ■' k =0 

The second term on the right tends to 0 when n goes to infinity, since n/ie(C/„) t and fc„ oo when 


n goes to OO. For the first term we assume first that (1 — e) ^ a. Then, 

/ dn i \ dn 


i—l k—Q 




• £)a 


-lli+l 


2 = 1 
dn 


(1 — <s) — 


2 = 1 


Since, as we have seen above, lim„^oo /?„ ~ * = 0, then lim„_>oo ^ — 0- Now, if 

(1 — e) 7 ^ Of, we have 


e/3n ( 1 + ^ ^(1 - efa j = em{Ako) fa^^" + dna‘^" + 0. 

V 7 = 1 i>—n / ^ / 


2=1 k—0 


Convergence of 

As in equation (l?51) we consider By (1^51) we have 

P(aW) 

P(H„) Me(C/„) ^ 

This means that the Extremal Index '§ is given by the noise level e. 


e =■. d. 
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